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Apparatus for Gas Filter Correlation Radiometry and Methods for 2- 
Dimensional and 3-Dimensional Atmospheric Sounding 



Related Cases 

This application claims the benefit of U.S. provisional application 60/255,225 
filed December 12, 2000, and U.S. provisional application 60/257,675 filed 
December 20, 2000 and which is hereby incorporated by reference into this 
application. 

Background of the Invention 

1. Field of the Invention 

This invention relates generally to instruments for measuring atmospheric or 
related gases by measuring infrared radiation. More particularly, this invention 
utilizes a plurality of gases that absorb or filter infrared radiation at frequencies (or 
equivalently at wavelengths) characteristic of each gas; the filtered radiation is 
passed to an infrared camera with associated electronics and signal processing, 
which compare the filtered radiation to unfiltered radiation to derive information about 
the atmospheric or related gases under study. 

2. Description of the Prior Art 

Gas filter correlation techniques have been used, in particular from 
spaceborne platforms, to perform atmospheric temperature sounding and to study 
atmospheric trace gases such as HCI, HF, CH 4 , and NO in the stratosphere. The 
gas correlation technique has also been employed from spaceborne platforms to 
study trace gases, especially CO, in the troposphere. Laser back scatter has been 
used to image distributions of CH 4 and related gases. Infrared cameras have been 
used extensively to produce images of infrared radiation sources. Fourier transform 
infrared (FTIR) instruments have been used to measure gas concentrations in the 
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atmosphere remotely. In situ air sampling instrumentation has also been used to 
measure gas concentrations in the atmosphere. 

While infrared cameras produce remotely sensed images of temperature 
distributions, they mostly measure surface temperatures, with limited information 
about air temperatures and even less information about gas concentrations. Air 
sampling instrumentation can measure gas concentrations accurately, but produces 
very limited information about spatial distributions. FTIR instrumentation can 
measure a large variety of chemical species simultaneously, but is rarely used in a 
true imaging system because of the enormously high data rates that result. Gas 
filter correlation techniques have been combined with imaging systems in the past, 
but generally with space based scanning imaging systems used only for low spatial 
resolution measurements. Laser back scatter systems require an active illumination 
source and are limited to wavelengths that can be produced efficiently by lasers, or 
by related coherent radiation sources. 

Pertinent prior art includes the HALOE, MAPS and MOPITT satellite 
instruments, and the MATR and MOPITT-A aircraft instruments. MAPS is described 
by H. G. Reichle, et al., in "Middle and Upper Tropospheric Carbon Monoxide Mixing 
Ratios as Measured by a Satellite-Borne Remote Sensor During November 1981", J. 
Geophys. Res., 91, pages 10,865-10,887 (Sept 20, 1986). MOPITT is described by 
J.R. Drummond, in "Measurements of Pollution in the Troposphere (MOPITT)", in 
The Use of EOS for Studies of Atmospheric Physics, pages 77-101, edited by J. C. 
Gille and G. Visconti, published by North-Holland (1992). MATR is described at 
http://www.eos.ucar.edu/matr/Welcome html. MOPITT-A is described at 
http://www.atmosp.physics.toronto.ca.MOPPITT .home.html. 

Additional prior art includes the following: J. X. Warner, et al. "Cloud detection 
and clearing for the Earth Observing System TERRA satellite Measurements of 
Pollution In The Troposphere (MOPITT) experiment", Appl. Opt., Vol. 40, No.8, 
March, 2001; Gas Correlation Spectroscopy, Opoto-Knowledge: The Source for 
Spectral Imaging - Press, http://www.techexpo. com/WWW/opto-knowledge/gas- 
corr.html, Oct. 3, 2000; J. Sandsten, et al., Gas Imaging Using Gas-Correlation 
Spectroscopy, Atomic Physics, LTH, Lena. Roos@fysik.lth.se,990603, http:// 
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www.atom-fysik.lth.se/AFDOC/ Progrep978/c3.htmi, Oct. 3, 2000; Mark Smith, 
MATR, MOPITT Airborne Test Radiometer, http://www.eos .edu .matr.Welcome.html, 
Sept. 20, 2000; Measurements of Pollution in the Troposphere, MOPITT, 
http.//www.atmosp .physics. utoronto.ca/MOPITT/home.html, Sept. 20, 2000; 
Measurement of Pollution in the Troposphere, MOPITT Project Home Page, 
http://www.eos .ucar.edu/mopitt/ home.html, (Atmospheric Chemistry Division, 
National Center for Atmospheric Research, Boulder, CO.), Sept. 20, 2000; 
Measurements of Pollution in the Troposphere, MOPITT Overview, http://www 
.atmosp.physics.utoronto.ca/MOPITT /overview, htm I, Sept. 20, 2000; Measurement 
of Pollution in the Troposphere (MOPITT), http://www.acd.ucar .edu/asr99/ 
MOPITT.html, Measurement of Pollution in the Troposphere (MOPITT) Program, 
Sept. 20, 2000; C. J. Clerbaux, et al., "Assimilation of carbon monoxide measured 
from satellite in a three-dimensional chemistry-transport model", in press, J. 
Geophys. Res., 2000; D. P. Edwards, et al., "Improvements to the correlated-k 
radiative transfer method: Application to satellite infrared sounding", J. Geophys. 
Res., Vol. 105, No. D14, pp. 18,135- 18,156, 2000; "ERAST: Scientific Applications 
and Technology Commercialization", Proc. of the NASA ERAST Exclusive Preview 
at NASA Dryden FRC, October 13, 1999, compiled by J. D. Hunley, et al., NASA/CP- 
2000-209031, Sept. 2000; B. V. Khattatov, et al., "On applications of photochemical 
models to the design of measurement strategies", submitted to Geophys. Res. Lett., 
2000; B. V. Khattatov, et al., "Assimilation of satellite observations in global 
chemistry-transport models", J. Geophys. Res., Vol. 105 , pp. 29135, 2000; C. D. 
Rodgers, "Inverse Methods for Atmospheric Sounding," Chapt. 2, (F.W. Taylor 
(University of Oxford) World Scientific, 2000); G. L. Stephens, et al., "The 
Department of Energy's Atmospheric Radiation Measurement (ARM) Unmanned 
Aerospace Vehicle (UAV) Program", Bull. Of the Amer. Meteor. Soc, Vol. 81, No. 
12, pp. 2915-2937, Dec. 2000; George V. Bailak, et al., MOPITT Airborne Validation 
Instrument: MOPITT-A, SPIE, (University of Toronto, Department of Physics, 
Toronto, ON, Canada) (Institute for Space and Atmospheric Studies, University of 
Saskatchewan, Saskatoon, SK, Canada), (Part of the SPIE Conference on Optical 
Spectroscopic Techniques and Instrumentation for Atmospheric and Space 
Research III, Denver, CO., July 1999), Vol. 3756, pp. 468-474; M.W. Smith, et al., 
"Remote sensing of atmospheric carbon monoxide with the MOPITT Airborne Test 
Radiometer (MATR)" Proc. SPIE ,Vol. 3756, pp. 475-485, July, 1999; D. P. Edwards, 
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et al., "Radiative transfer modeling for the EOS Terra satellite Measurement of 
Pollution in the Troposphere (MOPITT) instrument," J. Geophys. Res., Vol. 104, No. 
D14, pp. 16,755-16,775, July 27, 1999; J. F. Lamarque, et al., "Assimilation of 
Measurement of Air Pollution From Space (MAPS) CO in a global three-dimensional 
model", J. Geophys. Res., Vol. 104, pp. 26,209-26,218, 1999; L. Pan, et al., 
"Retrieval of tropospheric carbon monoxide for the MOPITT experiment", J. 
Geophys. Res., Vol. 103, No. D24, pp. 32,277-32,290, Dec. 27, 1998; P. F. Levelt, 
et al., "Assimilation of MLS ozone measurements in the global three-dimensional 
chemistry-transport model ROSE", Geophys. Res. Lett., 25, pp. 4493-4496, 1998; 
Mark W. Smith, et al., The Measurements of Pollutants in the Troposphere 
(MOPITT) Airborne Test Radiometer (MATR), The Faith Observer (Atmospheric 
Chemistry Division, National Center for Atmospheric Research, Boulder, CO.), 
http://eospso.gsfc.nasa.gov/eos _observ/ 7_8_98/p21.html, July/August 1998, Vol. 
10, No. 4; B. T. Tolton et al., Characterization of the Length-Modulated Radiometer, 
Applied Optics, Vol. 36, No. 22, pp. 5409-5420, August 1, 1997; M. W. Smith, 
"Method and results for optimizing the MOPITT methane bandpass", Appl. Opt., Vol. 
36, No. 98, pp. 4285-4291, June 20, 1997; Liwen Pan, et al., Satellite Remote 
Sensing of Tropospheric CO and CH 4 : Forward Model Studies of the MOPITT 
Instrument, Applied Optics, Vol. 34, No. 30, pp. 6976-6988, October 20, 1995; E. 
Andersson, et al., "Use of cloud-cleared radiances in three/four-dimensional 
variational data assimilation," Q. J. R. Meteorol. Soc, 120, pp. 627-653, 1994; J. M. 
Russell III, et al., "The Halogen Occultation Experiment", J. Geophys. Res., Vol. 98, 
No. D6, pp.1 0,777-10,797, June 20,1993; C. Grassotti, et al., "A Study of Satellite 
Emission Computed Tomography", Advances in Remote Sensing Retrieval Methods, 
edited by A. Deepak, H. E. Fleming, and J. S. Theon, A. Deepak Publishing, pp. 23- 
34, 1989; J. R. Drummond, Novel Correlation Radiometer: The Length Modulated 
Radiometer, Applied Optics, (Colorado School of Mines, Physics Department, 
Golden, CO) Vol. 28, No. 13, pp. 2451, July 1, 1989; Henry G. Reichle, Jr., et al., 
Feasibility of Determining the Vertical Profile of Carbon Monoxide from a Space 
Platform, Applied Optics, Vol. 28, No. 11, pp. 2104-2110, June 1, 1989; U.S. Patent 
No. 4,641,973 (issued Feb. 10, 1987) and U.S. Patent No. 4,4780,613 (issued Oct. 
25 1988); C. L. Wyatt, Radiometric System Design, MacMillan Publishing, Chap. 8, 
pp. 109-113, 1987; A. C. Lorenc, "Analysis methods for numerical weather 
prediction", Q. J. R. Met. Soc, 112, pp. 1177, 1986; H. E. Fleming, "Temperature 
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Retrievals Via Satellite Tomography", Advances in Remote Sensing Retrieval 
Methods, edited by A. Deepak, H. E. Fleming, and M. Chahine, A. Deepak 
Publishing, pp. 55-69, 1985; H. E. Fleming, "Satellite Remote Sensing by the 
Technique of Computed Tomography", J. Appl. Meteor., Vol. 21, pp. 1538-1549, Oct. 
1982; W. L. Smith, et al. "The Use of Interferometric Radiance Measurements for 
Sounding the Atmosphere", J. of Atmos. Sci., Vol. 36, pp. 566-575, April 1979; C. 
B. Ludwig, et al, Measurement of Air Pollutants from Satellites. 1: Feasibility 
Considerations, Applied Optics, Vol. 13, No. 6, pp. 1494-1509, June 1974; D. E. 
Burch, et al., Instrument to Monitor CH 4 , CO, and C0 2 in Auto Exhaust, Philco-Ford 
Corp, prepared for the U.S. Environmental Protection Agency, Oct. 1973; P. Ellis, et 
al., Remote Sounding of Atmospheric Temperature from Satellites, IV. The Selective 
Chopper Radiometer for Nimbus D, Proc. Roy. Soc. Lond. A. Vol. 334, pp. 149-170 
(1973); F. W. Taylor, et al., Radiometer for Remote Sounding of the Upper 
Atmosphere, Applied Optics, Vol 11, No. 1, pp. 135-141, Jan. 1972; J. T. Houghton, 
et al., Remote Sounding of Atmospheric Temperature from Satellites, I. Introduction, 
Proc. Roy. Soc. Lond. A., (Clarendon Laboratory, University of Oxford, J. J. 
Thomson Physics Laboratory, University of Reading), (Communicated by P.A. 
Sheppard, F.R.S., April 14, 1970), Vol. 320, pp. 23-33 (1970); P.G. Abel, Remote 
Sounding of Atmospheric Temperature from Satellites, II. The Selective Chopper 
Radiometer for Nimbus D, Proc. Roy. Soc. Lond. A. Vol. 320, pp. 35-55 (1970); D. Q. 
Wark, et al., "Indirect Measurements of Atmospheric Temperature Profiles from 
Satellites:", Mon. Weath. Rev., Vol. 94, No. 6, pp. 351-362, June 1966; L. D. 
Kaplan, "Inference of Atmospheric Structure from Remote Radiation Measurements", 
J. Opt. Soc. Am., 49, pp. 1004-1007, 1959; J. I. F. King, "The Radiative Heat 
Transfer of Planet Earth", Scientific Uses of Earth Satellites, edited by J. A. van 
Allen, The University of Michigan Press, pp. 133-1 36, 1958; Atmospheric Absorption, 
Field Measurements of Atmospheric Transmittance, 5-88, Fig. 5-31 Atmospheric 
Transmission at Sea Level over a 0.3 km path; L. Pan, et al., Analysis and 
Characterization of the Retrieval Algorithm for Measuring Tropospheric Co Using the 
Mopitt Instrument, SPIE (National Center for Atmospheric Research, Boulder, CO.), 
Vol. 2830/159; F. W. Taylor, Chapter 3, Pressure Modulator Radiometry, 
Spectrometric Techniques, (Department of Atmospheric Physics, Oxford University, 
Oxford, England),Vol. Ill; Boyd T. Tolton, et al., Calibration of a Length Modulated 
Radiometer, SPIE, (University of Toronto, Department of Physics, Toronto, Ontario, 
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Canada),Vol. 2830, pp. 253-263; J. M. Russell III, et al, Global monitoring of 
stratospheric halogen compounds from a satellite using gas filter spectroscopy in the 
solar occupation mode Appl. Opt., Vol. 15, No. 3, pp. 607-612, March 1977; U. 
Goers, et al., "A PPLN-OPO-backed backscatter absorption gas imaging (BAGI) 
system and its application to the visualization of fugitive gas emissions," SPIE 
(Sandia National Laboratories, Livermore, CA), (Department of Physics, University of 
Dayton, Dayton, OH), (Laser Imaging Systems, Inc., Punta Gorda, FL), Vol. 3758, 
pp. 172-179; K. W. Aniolek, et al., "Trace gas detection in the mid-IR with a compact 
PPLN-based cavity ringdown spectrometer," SPIE (Sandia National Laboratories, 
Livermore, CA and Albuquerque, NM), (Department of Physics, University of Dayton, 
Dayton, OH) Vol. 3758, pp. 62-73. 



The subject invention is a gas filter correlation camera and related methods. 
It is based on a novel combination of two techniques: the use of gas filter correlation 
radiometry in combination with an infrared camera. The camera comprises a 
rotatable filter wheel and an infrared camera. The filter wheel is provided with a 
plurality of chambers, each being filled with a different gas. Each gas corresponds to 
a gas of interest in the gas volume under study, usually the atmosphere. The 
camera can be used to produce three dimensional soundings of gas distributions by 
combining multi-angle viewing with tomographic, or related, reconstruction and 
retrieval techniques. 

Infrared light emitted by the gas volume under study is passed through each 
of the chambers as the filter wheel rotates, selectively filtering the light in a manner 
consistent with the infrared radiation absorption characteristics of the various gases. 
By correlating the position of the filter wheel with the receipt by the infrared camera 
(and related data-processing equipment) of filtered infrared radiation from the filter 
wheel, and by analyzing that light for each gas in the filter wheel, the temperature, 
concentration and other information for each of the corresponding gases in the gas 
volume under study may be obtained. 



Summary of the Invention 
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One unique aspect of this invention is to combine high spatial resolution 
images with high specificity to a suite of gases. In principle, any gas that is relatively 
stable and that has strong infrared absorption features could be imaged. One suite 
of gases that would be of strong commercial interest is comprised of gases of 
hydrocarbons, including substituted and unsubstituted saturated and unsaturated 
hydrocarbons, and particularly including straight-chain compounds like CH 4 
(methane) and C 3 H 8 (propane). Another important suite of gases is comprised of 
CO, C0 2 , and possibly 0 3> all of which are important by products of combustion. 
Another important suite of gases includes species that are emitted by volcanoes, 
such as HCI, HF, and S0 2 . Surface temperature and air temperature (possibly using 
N 2 0 or possibly using a second C0 2 band) would be measured in addition to the 
gases listed. Several of the gases listed, in particular CH 4 , C0 2 , and N 2 0, are 
important "greenhouse" gases. 3-D volume soundings could be created by means of 
tomographic reconstruction techniques if the instrument is flown over a given scene 
using different look angles. 

A principal objective of the invention is to produce two dimensional images 
(and/or three dimensional soundings) of certain geophysical parameters. These 
include surface temperature, air temperature, and the various suites of gases listed 
above. 

Brief Description of the Drawings 

FIG. 1 is a schematic diagram of a first gas correlation camera of the 
invention. 

FIG. 2 is a schematic diagram of a second gas correlation camera of the 
invention. 

FIG. 3 is a set of transmittance plots for gases used in the cells of a gas 
correlation camera of the invention. 
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FIG. 4 is a depiction of a first embodiment of a gas correlation filter wheel with 
multiple gas cells used to filter infrared spectra received by a gas correlation camera 
of the invention. 

FIG. 5 is a depiction of a second embodiment of a gas correlation filter wheel 
with multiple gas cells used to filter infrared spectra received by a gas correlation 
camera of the invention. 

FIG. 6 is a depiction of a third embodiment of a gas correlation filter wheel 
with multiple gas cells used to filter infrared spectra received by a gas correlation 
camera of the invention. 

FIG. 7 is a depiction of a fourth embodiment of a gas correlation filter wheel 
with multiple gas cells used to filter infrared spectra received by a gas correlation 
camera of the invention. 

FIG. 8 is a depiction of the gas correlation filter wheel of FIG. 6 and FIG. 9 
illustrates a method of rotating the filter wheel so as to maximize the time that each 
gas cell is receiving infrared light for passage to the camera receptor. 

FIG. 9 is a schematic diagram of a third gas correlation camera of the 
invention. 

FIG. 10 is a two-dimensional cross sectional view illustrating multi-angle 
viewing geometry for tomographic sounding using a gas correlation camera of the 
invention. 

FIG. 1 1 is an example of the three-dimensional density distribution derived 
using the tomographic sounding methodology of FIG. 10. 

Detailed Description of Preferred Embodiments 

This invention, when constructed in a physical embodiment, would have 
many potential uses. Four of the most important would be detecting emissions of 
methane (i.e. natural gas) and related hydrocarbons, monitoring and studying 
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wildfires, monitoring and studying air pollution (including greenhouse gases), and 
monitoring and studying volcanic emissions. 

The most likely commercialization of the invention would be for use in 
detecting emissions of natural gas and other hydrocarbons. The ability to conduct 
aerial surveys and produce high spatial resolution maps showing the distribution of 
natural gas (or methane) in the atmosphere above production fields and around 
storage and transportation facilities would be invaluable in tracking leaks or 
unintentional releases of this key energy source. Extending the detection capability 
to related volatile hydrocarbons, such as propane, would open the doors for still 
more useful applications. 

Another likely commercialization of the invention would be for use in 
monitoring wildfires and providing support for fire fighting efforts. For example, the 
invention could be used to provide early warning of the location of small, smoldering 
fires, thus allowing these fires to be quenched before they become full blown 
problems. For this purpose, the invention could be used to make aerial surveys over 
high risk areas, using elevated levels of carbon monoxide as a signature for a 
smoldering fire. This application would require full technical development of the 
instrument, but probably little development of scientific theory. The invention would 
probably not even have to provide fully calibrated (or quantitatively accurate) results, 
but rather provide ready discrimination of levels that are elevated compared to the 
normal background. Because the instrument uses thermal infrared radiation, it could 
be used day and night. 

For another example, the invention could be used to provide predictive 
support concerning the probable behavior of wild fires. This application would 
require extensive development of the underlying scientific theories that predict the 
behavior of fires, and would most likely require quantitatively accurate results, in 
order for these results to be incorporated in a useful way in predictive computer 
models. The underlying assumption is that the quantities measured (namely surface 
temperature, air temperature, CO and C0 2 concentrations) are strongly related to 
combustion efficiency, and therefore could be useful predictors of the evolution of a 
fire. The purpose would be to provide information that could be used to direct the 
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efforts of fire fighters for greatest efficiency, and also to provide early warning of 
possible rapid growth of a fire. 

A third possible role for the invention would be for monitoring and studying air 
pollution, conducted on either a regional basis using an airborne version of the 
instrument, or on a global basis using a spaceborne instrument. Both CO and C0 2 
are primary combustion by-products produced by industrial sources, by 
transportation sources, and by fires, both naturally occurring and started by humans. 
CO, being chemically reactive in the atmosphere, has a great impact on tropospheric 
oxidation cycles. C0 2 , as is well known, plays an important role in moderating the 



thermal environment of the earth and it's atmosphere, and hence has a big impact 
fiy on climate. By using N 2 0, or a second C0 2 band, for temperature sounding, this 
f j instrument could potentially provide a relatively simple means of making global 
W measurements of C0 2 concentrations. 

m 

O A fourth role for the invention would be to monitor and study volcanoes. 

nfi 

□ Operating from a fixed site, the instrument could ultimately provide simultaneous 

C| images (in real time) of air temperature and concentrations of HCI, HF and S0 2 in 

H the plume that issues from a volcanic vent. This could be useful information for 
predicting imminent eruptions. 

The equations that describe radiation transfer in the atmosphere are the basis 
for calculating atmospheric radiance and its sensitivity to changes in trace gas and 
temperature fields. These equations can be cast into several different forms. 
Equation 1 is one straightforward version: 



Zobs 



f u (v,zobs) = I v (u,zs)T(o,zs,zobs)+ J B(v,T(z)) Hki(p,z)Pi(z) T(V,Z,Zobs)dz 



(1) 



2s 



where 

v = wavenumber, with units of cm 1 

/ v (v 5 z obs ) = spectral radiance at v at observer altitude z obs , with units of W/(m 2 .sr~ ! .cm 1 ) 
7 v (v, z s ) = spectral radiance leaving the surface, also with units of W/(m 2 \sf k \cm l ) 
z(v, Zs, z obs ) = transmittance at v between z s and z obs , which is a dimensionless quantity 
B(y, T(zJ) = Planck function for temperature T at altitude z 5 with units of W/(m A \sf k \cm l ) 
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k,{v 9 z) = monochromatic absorption coefficient of i* absorber, units l/(molecules.cm 2 ) 
pi (z) = number density of i* absorber at altitude z, with units of molecules, cm' 3 
r(v, z, z 0 t,s) = transmittance at wavenumber v from altitude z to the observer at z Q bs 

The radiance that is detected by the n lh instrument channel is the product of 
atmospheric spectral radiance and instrument spectral response for that channel, 
integrated over wavenumber, as shown by Equation 2: 

In 9 obs = il u (o,zobs)R n (v)dv 

(2) L> = — OO 

where 

/ obs = band-integrated radiance observed by n 1 * 1 channel, with units of W/(m 2 .sr l ) 
I v (y, zobs) = spectral radiance at wavenumber v at the observer altitude z obs (from Eq. 1) 
R n (v) = spectral response function of n th instrument channel at wavenumber v 

The remote sensing process measures electromagnetic radiation, and the 
profiles of either temperature or trace gases are estimated, or retrieved, based on 
the amount of electromagnetic radiation that is measured at various wavelengths 
(multi-spectral approach) or from various angles (multi-angle approach). 

The equations that predict radiometric performance for a specific instrument 
can also be cast into several different forms. The following Equation 3 has been 
slightly rearranged from forms presented by Wyatt (Wyatt, 1987): 

SNR _ 2i^fMy^ zoW)dv 

P u-~co 

where 

SNR = signal-to-noise ratio, which is a dimensionless quantity 
D* = spectral detectivity for a specific detector, with units of cm 
v = wavenumber, with units of cm 1 

x e = efficiency of optics, which is a dimensionless quantity 

A d = area of detector (for a detector which is also the system field stop), with units of cm 2 
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/= frequency bandpass for electronic noise, with units of Hz 
F = relative aperture, which is a dimensionless quantity 
P = chopping factor, which is a dimensionless quantity 

The remaining terms in the integral are identical to those of Eq. 2. Eqs. 1 and 
2 are used to calculate the radiance that is expected for a given scene. Eq. 3 is 
used to predict the signal-to-noise ratio that is expected for a specific instrument that 
observes this scene. 

A representative example of a gas correlation camera of the invention is 
shown in FIG. 1. Infrared radiation 101 enters respective cells of filter wheel 102. 
The gas in each cell is different. The gas in each cell absorbs in a characteristic 
manner a particular band of the infrared radiation spectrum, except for at least one 
cell that does not absorb infrared radiation. Transmittance plots for gases used in 
the cells of a gas correlation camera of the invention are set forth in FIG. 3. 

Imaging optics 103 produce a desired field of view and correct optical 
misalignments of the infrared radiation samples as they exit the filter wheel cells. 
Infrared radiation detector 104 is a 2-D array detector with particular sensitivity to the 
infrared radiation band under study. Detector 104 receives and detects the infrared 
radiation emanating from the imaging optics 103. Signal conditioning electronics 105 
condition the electrical output of detector 104 before it is received by data acquisition 
computer 106. The computer compares the infrared radiation sample received from 
the non-absorbing wheel cell to the infrared radiation samples received from the 
other wheel cells to determine the spectral profile and energy content of the infrared 
radiation absorbed by each of the respective gases in the cell. This information 
permits determination of physical characteristics of the gas volume under study. 
Data generated by data acquisition computer 106 is displayed on real-time display 
monitor 107. 

The rotating gas wheel 1 02 is a key element of the correlation camera. 
Preferred embodiments include two distinct implementations of rotating gas wheel 
102. The wheel may be located either in front or inside the optics of infrared camera 
104. The advantage of locating wheel 102 in front of infrared camera 104 is 
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simplicity, while the main disadvantage is the fact that wheel 102 must then have its 
maximum size, as the first lens in a camera is generally the largest optical aperture. 
Locating the wheel in the midst of the camera optics would be more complicated, but 
could allow a smaller wheel to be used. This approach is discussed below in 
connection with FIG. 9. 

Another embodiment of a gas correlation camera of the invention is shown in 
FIG. 2. Filter wheel 202 of this embodiment comprises some chambers containing 
gases that absorb at one distinct band of the infrared radiation spectrum and other 
chambers containing gases that absorb at another distinct band of the infrared 
radiation spectrum, as well as a nonabsorptive cell. A first band could be in a range 
from about 2000 cm" 1 to about 2300 cm' 1 and a second band could be in a range 
from about 1300 cm' 1 to about 1360 cm" 1 . 

Infrared radiation 201 enters respective cells of filter wheel 202. Imaging 
optics 203 produce a desired field of view and correct optical misalignments of the 
infrared radiation samples as they exit the filter wheel cells as the filter wheel 202 
rotates to expose each successive cell to the imaging optics 203. Beam splitter 204 
splits the infrared radiation data from imaging optics 203 into two beams. One beam 
enters infrared radiation detector 205, which is an infrared camera with particular 
sensitivity to the distinct infrared radiation band absorbed by one set of the gases in 
filter wheel 202. The other beam from beam splitter 204 enters infrared radiation 
detector 209, which is an infrared camera with particular sensitivity to the distinct 
infrared radiation band absorbed by the other set of the gases in filter wheel 202. 
The other components of the system of FIG. 2 performs the same functions in the 
same way as the corresponding components of FIG. 1 . 

Infrared radiation detectors 205 and 209 are 2-D array detectors with 
particular sensitivity to the infrared radiation band under study. They receive and 
detect the infrared radiation emanating from the imaging optics 203 and beam 
splitter 204. Signal conditioning electronics 206 and 210 condition the electrical 
outputs of detectors 205 and 209 before they are received by data acquisition 
computers 207 and 21 1 . Each computer compares the infrared radiation samples 
received from the non-absorbing wheel cells to the infrared radiation samples 
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received from the other wheel cells to determine the spectral profile of the infrared 
radiation absorbed by each of the respective gases in the cell. This information 
permits determination of the physical characteristics of the gas volume under study. 
Data generated by data acquisition computers 207 and 21 1 is displayed on real-time 
display monitors 208 and 212. 

A first implementation of a rotating gas wheel is illustrated in FIG. 4. It is a 5- 
cell wheel 400 that is designed for use over the frequency range from about 2000 
cm" 1 to about 2300 cm" 1 . The body 401 of the cell pivots about axis 402 and 
contains 5 separate cells: 403, 404, 405, 406 and 407. The gases in each of these 
cells are as follows: 



1. 


Cell 403 = 


no gas (or optically transparent gas such as N 2 ); 


2. 


Cell 404 = 


N z O; 


3. 


Cell 405 = 


CO; 


4. 


Cell 406 = 


C0 2 ; and 


5. 


Cell 407 = 


band pass interference filter centered on O3 absorption features 



near 2100 cm' 1 ; note that Cell 407 actually does not contain Q 3 (very 



reactive) but rather a narrow band pass filter in lieu of O3. 

The gas cells will be separated by opaque, black sectors 408, so that a light 
interrupter or on/off chopper is an integral part of wheel 400. This wheel is useful for 
general studies of combustion sources. FIG. 3 illustrates the spectral absorption 
characteristics of these gases. The information that will be derived from radiometric 
measurements using each of these cells is as follows: 



1. 


Cell 403 = 


surface temperature; 


2. 


Cell 404 = 


air temperature; 


3. 


Cell 405 = 


CO amounts; 


4. 


Cell 406 = 


C0 2 amounts; 


5. 


Cell 407 = 


0 3 amounts. 



A second implementation of a rotating gas cell wheel of the invention is 
illustrated in FIG. 5. It is a 4-cell wheel 500 for use over a range from about 1300 
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cm' 1 to about 1360 cm' 1 . The body 501 of the wheel pivots about axis 502 and 
contains 4 separate cells: 503, 504, 505, and 506. The gases that will be contained 
in each of these cells are: 

1 . Cell 503 = no gas (or optically transparent gas such as N 2 ); 

2. Cell 504 = CH 4 ; 

3. Cell 505 = C0 2 ; 

4. Cell 506 = S0 2 . 

The information that will be derived from radiometric measurements using 
each of these cells is: 

1 . Cell 503 = surface temperature; 

2. Cell 504 = CH 4 amounts or air temperature; 

3. Cell 505 = C0 2 amounts or air temperature; and 

4. Cell 506 = S0 2 amounts. 



Again, the gas cells will be separated by opaque, black sectors 507, so that 
an on/off chopper is an integral part of wheel 500. This wheel will be useful for 
studies of strong sources of S0 2 in volcanoes. It could also be useful for studying 
sources of CH 4 or C0 2 . Which of these gases is used to sound for air temperature 
and which is treated as a species whose concentration is to be measured will 
depend on the particular application. 

A third implementation of a rotating gas cell wheel of the invention is a 5-cell 
design illustrated in FIG. 6. This wheel 600 has trapezoidal gas cells (603, 604, 605, 
606 and 607) rather than the round cells of the wheels of FIGS. 4 and 5. It rotates 
about pivot 602. The cells may be triangular rather than trapezoidal. 

The purpose of this design is to maximize the surface area of the cells on the 
wheel. This reduces noise and results in cleaner and more reproducible data in 
analysis of the signals generated by the CCD detector of the infrared camera. 
Again, the gas cells will be separated by opaque, black sectors 608, so that an on/off 
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chopper is an integral part of wheel 600, although the width of the chopper segments 
is maximized as much as possible. 

The wheel 600 of FIG. 6 can be configured similarly to the wheel 400 of FIG. 
4 so as to be useful over the frequency range from about 2000 cm" 1 to about 2300 
cm" 1 . The body 601 of the cell pivots about axis 602. The gases in each of the 
separate 5 are as follows: 



1. 


Cell 603 


= no gas (or optically transparent gas such as N2); 


2. 


Cell 604 


= N 2 0; 


3. 


Cell 605 


= CO; 


4. 


Cell 606 


= C0 2 ; and 


5. 


Cell 607 


= band pass interference filter centered on O3 absorption features 



near 2100 cm' 1 . 



This wheel is useful for general studies of combustion sources. The 
information that will be derived from radiometric measurements using each of these 
cells is as follows: 



1. 


Cell 603 


= surface temperature; 


2. 


Cell 604 


= air temperature; 


3. 


Cell 605 


= CO amounts; 


4. 


Cell 606 


= C0 2 amounts; 


5. 


Cell 607 


= O3 amounts. 



A fourth implementation of a rotating gas cell wheel of the invention is a 5-cell 
wheel 700 illustrated in FIG. 7. Here the surface area of the gas cells (703, 704, 
705, 706 and 707) has been enlarged to the maximum extent possible without 
compromising the structural integrity of the wheel itself, which rotates about pivot 
702. The chopper spindles 701 have correspondingly been reduced to the maximum 
extent possible. This wheel is probably not useful for the types of measurements 
discussed above because the presence of some significant chopper structure is 
important to making those measure accurately. This wheel would, however, be 
useful for measuring solar infrared radiation. 
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FIG. 8 is a depiction of the gas correlation filter wheel 800 of FIG. 6 illustrating 
a method of rotating the 5-cell filter wheel so as to maximize the time that each gas 
cell ( 803, 804, 805, 806 and 807) is receiving infrared light for passage to the 
infrared camera. The purpose of this method is to increase the efficiency and 
throughput of the filter wheel and hence the gas correlation camera as whole. 

To this end, the wheel rotates, but not at a steady rate. The point is to reduce 
the time that the imaging optics see the chopper spindles 801 and correspondingly 
increase the time that the optics see each of the cells. Rotational velocity 81 1 
(shown as a vector between reference lines 809 and 810) is increased when a 
chopper spindle 801 is blocking the imaging optics. Rotational velocity 812 (shown 
as a vector between reference lines 808 and 809) decreases when a gas cell is seen 
by the imaging optics. The net effect is to minimize the time that the imaging 
electronics spend looking at chopper spindles and maximize the time spent looking 
at and deriving useful data from gas cells. 

FIG. 9 is a schematic diagram of a third gas correlation camera of the 
invention in which a rotating gas wheel is positioned inside the optics of an infrared 
camera. The advantage of locating the wheel in front of the infrared camera is 
simplicity, while the main disadvantage is the fact that the wheel must then have its 
maximum size, as the first lens in a camera is generally the largest optical aperture. 
Locating the wheel in the midst of the camera optics would be more complicated, but 
could allow a smaller wheel to be used. 

A gas correlation camera 900 locating the wheel in the midst of the camera 
imaging optics is illustrated in FIG. 9. Imaging optics 902 comprise a part of the 
overall imaging of the system. Imaging optics 907 complete the imaging optics of 
the system. Schematic 908 comprises IR detector 909 and conditioning electronics 
911. 

The key to this design is the location of relatively small gas filter wheel 903 
between imaging optics portion 902 and portion 907. As in designs previously 
discussed, the cells of the wheel are successively positioned in the optical axis of the 
imaging optics 902 and 907 so as successively to filter incoming raw infrared 
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radiation 901. After filtration, infrared radiation detected by detector 909 and a 
corresponding electrical signal 910 is passed to signal conditioning electronics 911. 

Wheel 903 is fitted with gear teeth on its exterior. These gear teeth 
correspond to gear teeth on drive gear 904, which is driven by motor 906 through 
shaft 906. There are many alternative drive mechanisms known to one of skill in the 
art. For example, drive gear 904 could instead be a pulley for driving wheel 903. 

Figure 10 illustrates a system for conducting tomographic sounding using one 
or more gas correlation cameras of the invention. In this embodiment, the 
correlation camera(s) are carried either by one or more conventional aircraft or by 
one or more unmanned aerospace vehicles ("UAVs") 1000 and 1001 and achieve 
overlapping multi-angle views of the atmosphere 1002 and 1003; these views form 
the basis for 3-D tomographic retrievals. This system will retrieve 3-D distributions of 
key tropospheric trace gases, as well as air and surface temperature. It also will 
take advantage of the unique capabilities of aircraft or UAVs as highly programmable 
and maneuverable observation platforms. While the system of FIG. 10 carries gas 
correlation cameras, other instrumentation certainly could be carried as part of the 
payload. Cloud and aerosol remote sensing capabilities would be especially 
valuable. 

While satellites are preferable to aircraft or UAVs for global observations, 
Aircraft and UAVs possess some unique advantages for high spatial and temporal 
resolution regional-scale and local-scale observations. An aircraft or UAV can fly in 
a holding pattern over a given scene to provide continuous temporal coverage of 
events for as long as the vehicle can remain airborne. The observation time can be 
extended indefinitely by using a small fleet of UAVs flying in relay fashion. Related 
to this is the possibility of long signal averaging periods to increase measurement 
precision. Aircraft and UAVs also offer tremendous flexibility in obtaining the multiple 
view angles that are required for tomographic sounding of the atmosphere. Finally, 
the aircraft or UAV flight altitude (up to tens of kilometers) will be much lower than a 
satellite orbital altitude (hundreds to thousands of kilometers). This will permit much 
higher spatial resolution in all dimensions. 
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There are two distinct factors that contribute to the vertical resolution that is 
attainable with multi-angle remote sounding of the atmosphere. This has been 
pointed out especially by other workers in the field. One factor is the change in 
shape, and upwards shift, of the weighting function as the view angle moves off 
nadir. This effect is actually similar to what occurs in the multi-spectral approach, as 
it is related to the rate of change of optical depth as a function of altitude. The 
second factor is the possibility of applying tomographic retrieval techniques to data 
that includes multiple views, acquired from multiple angles, of each portion of the 
atmosphere. For tomographic retrievals, the instrument's field of view, and not just 

2 the weighting functions, determines the vertical resolution. Tomographic sounding 
therefore has the potential to produce vertical resolution superior to that of 

it! conventional multi-spectral sounding, for which vertical resolution is limited, even for 

fSp the hypothetical case of perfectly monochromatic observations, to a couple of 

j^l kilometers at best and more typically to 8-12 km. 

S One embodiment of the system of FIG. 10 will use the general Bayesian 

D optimal estimation approach for reconstruction of the 3-D density distributions of 

atmospheric gases and temperatures from the multi-angle radiance observations. 
^ This reconstruction approach is not limited to a particular geometry, number, or type 

of observations, and it allows calculations of uncertainty of the derived density 

distributions, thus supplying an important statistical indicator of the quality of final 

results. 

In this technique crude a priori approximations of 3-D density distributions are 
corrected by observations in regions where observations are available and the 
calculated uncertainty of the derived distributions decreases accordingly in the 
corresponding locations. The more observations that are available the more 
accurate the resulting distributions become, and even in cases where very few 
observations are available the method still provides approximate estimates of spatial 
density distributions. The reliability of the estimates can be established from the 
supplementary uncertainty distributions. The approach can make use of previously 
(approximately) known properties of distributions of the tropospheric gases under 
consideration, such as the degree of spatial correlations in both the vertical and 
horizontal directions. These properties can be approximated from high-resolution 
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chemistry-transport model calculations or from past observations. Particular 
geometry, number, or type of observations does not limit the practical 
implementations of the approach. For instance, it is not necessary for the lines-of- 
sight of the observations shown in FIG. 10 to intersect in space as long as they 
sample roughly the same volume of interest. 

The optimal analysis equations can be presented as Equation 4: 



x a - resulting analyzed 3-D density distributions and temperatures arranged in a vector. 
x-a priori, approximate 3-D density distributions and temperatures arranged in a vector. 
y - observations of radiances arranged in a vector. 
K - Kalman gain matrix. 

B a - error covariance matrix for the analyzed (retrieved) 3-D distributions. 
B - error covariance matrix for the a priori distributions. 
O - observational error covariance. 

H- observational operator, i.e., conversion from gas densities and temperatures in the nodes 
of a regular grid to the observed radiances. This operator includes both the spatial 
interpolation and the non-linear forward radiation transfer model. 

H 9 - linearized observational operator, that includes linearized forward radiation transfer 
model. 



During the analysis the a priori distributions x obtained from chemistry- 
transport models or past observations are used to calculate radiances using the 
observational geometry and the forward model H. These calculated radiances, H(x), 
are then compared to the observed radiances y and spatial 3-D distributions x a of the 
target chemical species and temperatures are modified to minimize the difference 
between the calculated (H(x)) and the observed (y) radiances. The resulting field x a 
provides the best estimate of the true distribution of the target chemical species and 



x a =x +K(y-H(x)) 
K = B H ,T (HB H' T +0) 
B fl =(I-KH')B 



JSP, 



where 
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temperatures. Since operator H is non-linear, several iterations of this procedure 
might be required to reach convergence and the tangent linear model H' and the 
adjoint tangent linear model H' r of the utilized radiation transfer model need to be 
constructed. 

The off-diagonal elements of the background error covariance matrix B 
contain information about characteristic horizontal and vertical correlation scales of 
the tropospheric gas under consideration. They are usually computed as illustrated 
in Equation 5: 




where L h and L z are the characteristic correlation scales in the horizontal and vertical 
directions and Ar h and Ar z are, respectively, horizontal and vertical separations. 

A variation of the proposed technique has been applied successfully to the 
case of the assimilation of TOVS radiance in the ECMWF assimilation package. A 
similar approach has been applied to reconstructing spatial density distributions of 
various stratospheric and tropospheric trace gases using observations from satellite- 
based instruments . 

Fig. 1 1 is a three-dimensional plot of certain data collected by a NASA 
satellite. The satellite did not use the gas correlation camera of the invention. FIG. 
1 1 is included as a general illustration of a tomographic technique. 

Density distribution 1 101 is an example of a reconstructed isosurface of 
instantaneous three-dimensional density distribution of CO in the atmosphere. 
Density distribution 1 101 represents the density of CO gas in the atmosphere 
sampled in a window by the NASA satellite as it made three circles of the earth. 
Density distribution 1 101 is defined and scaled by Z-axis 1002, Y-axis 1003 and X- 
axis 1004. 
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Y-axis 1003 represents the latitudinal position of any given point of the density 
distribution. X-axis 1 104 represents the longitudinal position of any given point of the 
density distribution 1 004. Z-axis 1 002 represents the "model level," which 
corresponds to the position of any given point of the density distribution relative to 
the surface of the earth. Note that the density distribution has been flipped, so that 
the highest values on Z-axis 1002 represent points closest to the earth. The 
shadings on the density distribution plot 1 101 represent relative densities, with the 
highest density 1 105 being on the order of 200 ppb. 
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While the exemplary preferred embodiments of the present invention have 
been described in detail, those skilled in the art will recognize various changes, 
modifications, additions, and applications other than those specifically mentioned 
herein which fall within the spirit of this invention. 



Claims 

What is claimed is: 
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